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Abstract. With the Westerbork Synthesis Radio Telescope (WSRT), multi- frequency polarimetric images were 
taken of the diffuse radio synchrotron background in a ~ 5°x7° region centered on (l,b) = (161°, 16°) in the 
constellation of Auriga. The observations were done simultaneously in 5 frequency bands, from 341 MHz to 
375 MHz, and have a resolution of ~ 5.0' x 5.0' cosec <5. The polarized intensity P and polarization angle (j> show 
ubiquitous structure on arcminute and degree scales, with polarized brightness temperatures up to about 13 K. 
On the other hand, no structure at all is observed in total intensity I to an r.m.s. limit of 1.3 K, indicating 
that the structure in the polarized radiation must be due to Faraday rotation and depolarization mostly in the 
warm component of the nearby Galactic interstellar medium (ISM). Different depolarization processes create 
structure in polarized intensity P. Beam depolarization creates "depolarization canals" of one beam wide, while 
depth depolarization is thought to be responsible for creating most of the structure on scales larger than a beam 
width. Rotation measures (RM) can be reliably determined, and are in the range —17 RM 10 rad m -2 
with a non-zero average RMo w —3.4 rad m -2 . The distribution of RMs on the sky shows both abrupt changes 
on the scales of the beam and a gradient in the direction of positive Galactic longitude of ~ 1 rad m~ 2 per 
degree. The gradient and average RM are consistent with a regular magnetic field of ~ 1 ^iG which has a pitch 
angle of p = —14°. There are 13 extragalactic sources in the field for which RMs could be derived, and those 
have \RM\ ^ 13 rad m -2 , with an estimated intrinsic source contribution of ~ 3.6 rad m -2 . The RMs of the 
extragalactic sources show a gradient that is about 3 times larger than the gradient in the RMs of the diffuse 
emission and that is approximately in Galactic latitude. This difference is ascribed to a vastly different effective 
length of the line of sight. The RMs of the extragalactic sources also show a sign reversal which implies a reversal 
of the magnetic field across the region on scales larger than about ten degrees. The observations are interpreted in 
terms of a simple single-cell-size model of the warm ISM which contains gas and magnetic fields, with a polarized 
background. The observations are best fitted with a cell size of 10 to 20 pc and a ratio of random to regular 
magnetic fields B rarl / B reg w 0.7 ±0.5. The polarization horizon, beyond which most diffuse polarized emission is 
depolarized, is estimated to be at a distance of about 600 pc. 

Key words. Magnetic fields - Polarization - Techniques: polarimetric - ISM: magnetic fields - ISM: structure - 
Radio continuum: ISM 



1. Introduction 

The warm ionized gas and magnetic field in the Galactic 
disk Faraday-rotate and depolarize the linearly polar- 
ized component of the Galactic synchrotron emission. 
Structure of the warm ionized Interstellar Medium (ISM) 
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is imprinted in the polarization angle and polarized inten- 
sity by Faraday rotation and depolarization. Therefore, 
the structure in the polarized diffuse radio background 
provides unique information about the structure in the 
magneto-ionized component of the ISM. Using interferom- 
eters we can probe arcminute scales, which correspond to 
linear scales from a fraction of a parsec to several tens of 
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parsecs. Multi-frequency radio polarization observations 
of the Faraday-rotated synchrotron emission enable de- 
termination of the rotation measure (RM) distribution 
along many contiguous lines of sight. Therefore, this is a 
valuable method to estimate the non-uniform component 
of the Galactic magnetic field, weighted with electron den- 
sity. 

Although the magnetic field is only one of the fac- 
tors that shape the very complex multi-component ISM, 
it is generally thought to play a major role in the energy 
balance of the ISM and in setting up and maintaining the 
turbulent character of the medium. The non-uniform com- 
ponent of the magnetic field may influence heating of the 
ISM (Minter & Balser 1997), and provide global support 
of molecular clouds (e.g. Vazquez-Semadeni et al. 2000, 
Hcitsch et al. 2001). It also plays an important role in 
star formation processes (e.g. Shu 1985, Beck et al. 1996, 
Ferriere 2001). Furthermore, observations of the detailed 
structure of the magnetic field can provide constraints for 
Galactic dynamo models (Han et al. 1997). 

The strength and structure of the magnetic field in 
the warm ionized medium in the Galaxy can also be ob- 
tained from Faraday rotation measurements of pulsars, or 
polarized extragalactic sources. 

Among these, radio observations of pulsars take a spe- 
cial place because for them one can measure the RM and 
the dispersion measure (DM) along the same line of sight. 
The ratio of the two immediately yields the averaged com- 
ponent of the magnetic field along a particular line of 
sight (see e.g. Lyne & Smith 1989, Han et al. 1999). These 
studies, and analysis of the distribution of polarized ex- 
tragalactic sources (e.g. Simard-Normandin & Kronberg 
1980, Sofue & Fujimoto 1983), suggest that the regular 
Galactic magnetic field is directed along the spiral arms, 
with a pitch angle of p ~ 5 — 15°. There is evidence for re- 
versals of the uniform magnetic field inside the solar circle 
(e.g. Simard-Normandin & Kronberg 1979, Rand & Lyne 
1994), while there is still debate on reversals outside the 
solar circle (e.g. Vallee 1983, Brown & Taylor 2001, Han 
et al. 1999). Rand & Kulkarni (1989) used pulsar RMs 
to derive a value for the strength of the regular compo- 
nent of the magnetic field of B reg = 1.3 ±0.2 p,G. This 
result assumes a circularly symmetric large-scale Galactic 
magnetic field. The residuals with respect to the best-fit 
large-scale model were interpreted as due to a random 
magnetic field of 5 /iG. This number was derived by as- 
suming a scale length of the random field, modeled as a 
cell size, of 55 pc. Ohno & Shibata (1993) confirmed this 
result for the random magnetic field component for all cell 
sizes in the range of 10 - 100 pc. They used 182 pulsars in 
pairs and therefore did not have to make any assumptions 
about a large-scale Galactic magnetic field. 

Using extragalactic sources, most of which are double- 
lobed with lobe separations between 30" and 200", Clegg 
et al. (1992) found substantial fluctuations in RM on lin- 
ear scales of ~ 0.1 - 10 pc, which they could explain by 
electron density fluctuations alone. However, Minter & 
Spangler (1996) observed RM fluctuations in extragalac- 



tic source components that cannot be explained by only 
electron density fluctuations; they need an additional tur- 
bulent magnetic field of ~ 1 /iG to fit the observations. 
Observations of polarization of starlight (Jones et al. 1992) 
give much larger estimates of the cell size, of up to a kpc. 

Using pulsars and extragalactic radio sources for the 
determination of the RM of the Galactic ISM is clearly 
not ideal, because they only provide information in par- 
ticular directions which sample the ISM very sparsely. On 
the contrary, the diffuse Galactic radio background pro- 
vides essentially complete filling over large solid angles. 
Therefore, rotation measure maps of the diffuse emission 
can be produced that give information on the electron- 
density-weighted magnetic field over a large range of 
scales. The distribution of polarized intensity, when in- 
terpreted as mostly being due to depolarization, can yield 
estimates of several properties of the warm ISM such as 
correlation length, ratio of random over regular magnetic 
field and the distance out to which diffuse polarization can 
be observed. Early RM maps of the diffuse synchrotron 
emission in the Galaxy were constructed by Bingham & 
Shakeshaft (1967) and Brouw & Spoelstra (1976) who con- 
firmed a Galactic magnetic field in the Galactic plane. 
Junkes et al. (1987) presented a polarization survey of the 
Galactic plane at 4.9° < I < 76° and \b\ < 1.5° showing 
small-scale structure in diffuse polarization. The existence 
of polarization filaments at intermediate latitudes, with- 
out correlated structure in total intensity /, was discov- 
ered by Wieringa et al. (1993) at 325 MHz. Polarization 
surveys have been performed at frequencies from 1.4 GHz 
to 2.695 GHz (Duncan et al. 1997, Duncan et al. 1999, 
Uyamker et al. 1999, Landecker et al. 2001, Gaensler et 
al. 2001), mostly in the Galactic plane. 

In this paper, we discuss the results of multi-frequency 
observations at low frequencies around 350 MHz of a field 
in the constellation Auriga, in the second Galactic quad- 
rant (I = 161°, b = 16°). Due to the low frequencies, we 
probe low rotation measures that are predominant at in- 
termediate and high latitudes. This gives the opportunity 
to study the high-latitude RM without concrete objects 
such as HII regions or supernova remnants in the line of 
sight, and estimate structure in the ISM above the thin 
stellar disk. 

Section [5] contains details of the multi-frequency po- 
larization observations. In Sect. we analyze the obser- 
vations, in particular the small-scale structure in polar- 
ized intensity and polarization angle. Faraday rotation is 
discussed in Sect. where we also present the map of 
rotation measure. In Sect. [SJ depolarization mechanisms 
are described that cause the structure in P, and the con- 
straints that the observations provide for the parameters 
that describe the warm ISM. In Sect. H3 the polariza- 
tion properties of 13 polarized extragalactic point sources 
found in the Auriga region are discussed. In Sect. we 
discuss the information that our data provides on the 
strength and structure of the Galactic magnetic field. 
Finally, our conclusions are stated in Sect. [HI 
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2. The observations 

We used the Westerbork Synthesis Radio Telescope 
(WSRT) for multi-frequency polarimetric observa- 
tions of the Galactic radio background in a field 
in the constellation of Auriga centered on (a, S) 
(B1950) = (6 ft 20 m , 52°30 m ) (I = 161°, b = 16°). This 
field was observed in 8 frequency bands between 325 
and 390 MHz simultaneously, each with a band width of 
5 MHz. Due to radio interference and hardware problems 
only data in 5 of the 8 bands, viz. those centered on 341, 
349, 355, 360, and 375 MHz, could be used. 

The region in Auriga was observed in six 12hr periods, 
which resulted in a baseline increment of 12m. The short- 
est baseline obtained is 36m, and the longest is 2700m, 
which gives a maximum resolution of ~ 1 . A taper was 
applied to the (u, w)-data to increase the signal-to-noise 
ratio, so as to obtain a resolution of 5.0' x 5.0' cosec S = 
5.0' x 6.3' in all 5 frequency bands. As an interferometer 
has a finite shortest baseline, it is insensitive to large-scale 
structure. The shortest spacing of 36m in the WSRT con- 
stitutes effectively a high-pass filter for all scales above 
approximately a degree. 

The Auriga region was selected from diffuse polariza- 
tion maps that were produced as a by-product of the 
Westerbork Northern Sky Survey (WENSS, Rengelink 
et al. 1997), which is a single- frequency radio survey at 
325 MHz. The WENSS diffuse polarization maps contain 
several regions of high polarization which show conspic- 
uous small-scale structure in polarized intensity and po- 
larization angle. We reobserved two of those regions at 
multiple frequencies and with higher sensitivity to obtain 
rotation measure information. The first region in Auriga is 
described in this paper, and the other one in a forthcoming 
paper (Haverkorn et al. 2003a). 

The data reduction process is described in detail in 
Haverkorn (2002), and we only give a brief summary here. 
The observations were reduced using the NEWSTAR data 
reduction package. Polarized and unpolarized standard 
calibrator sources were used, where the absolute flux scale 
at 325 MHz is based on a value of 26.93 Jy for 3C286 
(Baars et al. 1977). From this value the flux scales of the 
other calibrator sources 3C48, 3C147, 3C345, and 3C303 
were derived. 

As the area to be mapped is larger than the primary 
beam of the WSRT, the mosaicking technique was used 
(Rengelink et al. 1997). In this observing mode, the ar- 
ray cycles through a preselected set of pointing positions 
a number of times during the 12 hour observation pe- 
riod, which reduces the instrumental polarization in be- 
tween pointing centers to below 1%. For the present re- 
gion we used 5x7 pointing positions separated by 1.25°. 
We only discuss the central ~ 7°x9° of the observed field 
of view, as the edges of the mosaic exhibit larger noise 
due to primary beam attenuation and increased instru- 
mental polarization. Maps were constructed of Stokes pa- 
rameters /, Q, U, and V. From Stokes Q and U, po- 
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Central position 


(Lb) = (161°, 16°) 


Size 


~ 7°x9° 


Pointings 


5x7 


Frequencies 


341, 349, 355, 360, 375 MHz 


Resolution 


5.0' x 5.0' cosec 5 = 5.0' x 6.3' 


Noise 


~ 4 mjy/beam (0.5 K) 


Conversion Jy - K 


1 mjy/beam = 0.127 K (at 350 MHz) 


Spacings 


Date Start time End time 
(UT) (UT) 


36m 


96/02/19 13:36 01:54 


48m 


95/12/24 17:28 05:41 


60m 


95/12/31 17:13 05:31 


72m 


96/01/07 16:37 04:55 


84m 


96/01/30 14:59 03:12 


96m 


96/01/09 16:33 04:46 



Table 1. Observational details of the Auriga field 



larized intensity P = y/ Q 2 + U 2 and polarization angle 
(f> — 0.5 arctan(J7/Q) 1 were derived. 

To avoid spurious solar emission coming in through po- 
larized side lobes, all observations were done at night, as 
shown in Table. ^ In addition, the observing period was 
close to solar minimum and in winter, minimizing iono- 
spheric Faraday rotation. The high polarized brightness 
that we observe (see Sect. 13.21) indicates that ionospheric 
Faraday rotation does not vary much during each of the 
12hr observations, because if it did, hardly any polarized 
intensity would have resulted. The differences between the 
average ionospheric Faraday rotation in the six 12hr obser- 
vations were determined from several point sources (in the 
Q- and [/-maps at ~ 1' resolution). The result is shown 
in Fig. ^ f° r 4 strong sources and 4 weaker sources. The 
strong sources have a signal-to-noise a = 2 - 3.5 for an 
individual 12hr observation, which corresponds to an un- 
certainty in polarization angle w 10°. The signal-to- 
noise of the weaker sources is a = 1.2 - 1.5 (o^ « 22°). 
From Fig. \I\we conclude that differential Faraday rotation 
between the six 12hr observations is negligible to within 
the uncertainties with which it can be measured from the 
point sources. 

We can estimate the RM contribution of the iono- 
sphere from the earth magnetic field and the electron den- 
sity in the ionosphere. The contribution of the ionospheric 
Faraday rotation is estimated as follows. The total elec- 
tron content (TEC) in the ionosphere above Westerbork 
at night, in a solar minimum, and in winter is minimal: 
TEC ~ 2.2 10 16 electrons m™ 2 (Campbell, private com- 
munication). The Auriga field at declination 5 — 52.5° is 
almost in the zenith at Westerbork at hour angle zero. 
If we assume a vertical component of the earth magnetic 
field of 0.5 G towards us, and a path length through the 
ionosphere of 300 km, the RM caused by the ionosphere 

1 Note that this notation, although much used, is actually in- 
correct for Q < 0. The correct notation is (f> — 0.5 angle (Q, U), 
or <f> = 0.5 arctan([//Q) for Q > and <j> = 0.5 arctan(C//Q) + 
0.5 * 7T for Q < 0. 



4 



M. Haverkorn et al.: Multi-frequency polarimetry around 350 MHz I 




ing the 2.7 K contribution of the cosmic microwave back- 
ground. Assuming a temperature spectral index of —2.7, 
the backgrounds at the lowest (341 MHz) and the high- 
est (375 MHz) frequency of observation are ~ 49 K and 
~ 38 K, respectively. Of these, approximately 25% is due 
to sources, as estimated from source counts (Bridle et al. 
1972), and assuming the spectral index of the Galactic 
background and the extragalactic sources to be identical. 
We thus estimate the temperatures of the diffuse Galactic 
background at 341 and 375 MHz to be 37 and 29 K, re- 
spectively. 



1 2 3 4 5 6 

observing night 

Fig. 1. Polarization angle variation over six 12hr observa- 
tions, for 8 extragalactic point sources. Four sources de- 
noted by the solid lines arc relatively strong (a — 2 - 3.5 for 
one 12hr measurement), the dotted lines are sources with 
a = 1.2 - 1.5. Little angle variation over nights points to 
constant ionospheric Faraday rotation. 

is 0.3 rad m -2 at hour angle zero. At larger hour angles, 
this is even less. So we expect the rotation measure val- 
ues given in this paper not to be affected by ionospheric 
Faraday rotation by more than 0.5 rad m~ 2 . 

2.1. Missing large scale structure 

An interferometer is increasingly insensitive to structure 
on larger angular scales due to missing small spacings. 
This means for the WSRT that scales above approxi- 
mately a degree cannot be detected. The Q and U maps 
are constructed so that in each mosaic pointing, the map 
integral of Q and U are zero. This leads to missing large- 
scale components in Q and U , and therefore erroneous de- 
terminations of P, (j> and RM . However, the distribution 
of RM is broad enough within one pointing. Therefore, 
the variation in polarization angle is so large that the av- 
erage Q and U are close to zero, so that missing large-scale 
components are negligible (Sect.JSJ. 

3. Total intensity and linear polarization maps 

3.1. Total intensity 

In the left-hand panel of Fig. [21 we show the map of to- 
tal intensity / at 349 MHz with a resolution of 5.0' x 6.3'. 
Point sources were removed down to the confusion limit 
of 5 mjy/beam. / has a Gaussian distribution around 
zero with width 07 = 1.3 K (while the noise, computed 
from V maps, is cry = 0.5 K), and is distributed around 
zero due to the insensitivity of the interferometer to scales 
> 1°. From the continuum single-dish survey at 408 MHz 
(Haslam et al. 1981, 1982), the /-background at 408 MHz 
in this region of the sky is estimated to be ~ 33 K 
with a temperature uncertainty of ~ 10% and includ- 



3.2. Polarized intensity 

In the right-hand panel of Fig. 2 we show the map of polar- 
ized intensity at 349 MHz, with a resolution of 5.0' x 6.3'. 
The noise in this map is ~ 0.5 K, and the polarized bright- 
ness temperature goes up to ~ 13 K, with an average value 
of- 2.3 K. 

The map shows cloudy structure of a degree to a few 
degrees in extent, with linear features that are sometimes 
parallel to the Galactic plane. In addition, a pattern of 
black narrow wiggly canals is visible (see e.g. the canal 
around (a, S) = (92.7°, 49° - 51°)). These canals are all 
one synthesized beam wide and are most likely due to 
beam depolarization. They separate regions of fairly con- 
stant polarization angle between which the difference in 
polarization angle is approximately 90° (±71 180°, n = 
1,2,3...) (Haverkorn et al. 2000). A change of 90° (or 
270°, 450° etc.) in polarization angle within one beam can- 
cels the polarized intensity in the beam; therefore these 
canals appear black in Fig. [21 Hence, the canals are not 
physical structures, but the reflection of specific features 
in the distribution of polarization angle. Such canals have 
also been detected by others at higher frequencies (e.g. 
Uyamker et al. 1999, Gaensler et al. 2001). It must be 
appreciated that the canals represent an extreme case of 
beam depolarization. Beam depolarization is much more 
wide-spread if the polarization angle varies on the scale 
of the beam or on smaller scales, but in general its effect 
is (much) less than in the canals. Therefore, with the ex- 
ception of the canals, beam depolarization does not leave 
easily visible traces in the polarized intensity distribution. 

An alternative explanation for depolarized regions ac- 
companied by a change in polarization angle of 90°is so- 
called depth depolarization or differential Faraday rota- 
tion. If the magnetic field in the medium is uniform, depo- 
larization along the line of sight causes depolarized regions 
at certain RM values (Sokoloff et al. 1998). However, in 
this scenario, it is hard to explain why the canals are all ex- 
actly one beam wide. Furthermore, the canals should shift 
in position with frequency, which is not observed. For an 
extended discussion on the possible causes of polarization 
canals, see Haverkorn et al. 2003b. 

There is no obvious structure in Stokes I, and if there 
is any, it does not appear to be correlated with the struc- 
ture in P. This is not typical for the Auriga field, but 
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96 94 92 90 88 96 94 92 90 88 

a (1950) (°) a (1950) (°) 

Fig. 2. Map of the total intensity 7 (left) and polarized intensity P (right) at 349 MHz, on the same brightness scale. 
The maps are smoothed to a resolution of 5.0' x 6.3'. White denotes high intensity, up to a maximum of ~ 13 K. 
Superimposed in the right map are lines of constant Galactic latitude, at latitudes of 13°, 16° and 19°. 



it is true in all fields observed so far with the WSRT at 
~ 350 MHz (Katgert & de Bruyn 1999, Haverkorn et al. 
2003a, Schnitzeler et al. in prep). 

We have estimated the amount of correlation between 
total and polarized intensity by deriving the correlation 
coefficient C, where C(f,g) of the observables f{x,y) and 
g{x,y) is defined as 



C(f,g) 



Ei,i(/ij -/) (m -g) 



T,i,j(fij ~ f ? T,iA9ij - a) 



(i) 



1 , N y , and N x and 



where / is the average value of f(x,y), /y 
summation is over i = 1, N x and j 
N y is the number of data points in the x- and y-directions, 
respectively. The correlation coefficients between P and 7, 
and between P and between 7 at different frequencies are 
given in Tabled Both the observed correlation coefficients 
at the maximum resolution of about 1', as well as those of 
the 5' resolution data are given. 

The high correlations between 7 in different frequency 
bands are expected if the visible structure in 7 is mostly 
due to faint sources. The C(I,I) < 1 is due to uncor- 
rected noise. However, the correlation between P and 7 
is very low, so the structure in polarization is not due 
to small-scale variations in synchrotron emission. Instead, 



Table 2. Correlation coefficients for correlations between 
P and 7, between 7 in different frequency bands and be- 
tween P in different frequency bands, for tapered data and 
full resolution data. 



freq band k 


res 


C(I k ,P k ) 


C{Pux,P k ) 


C(lMl,Ik) 


(MHz) 


(') 








341 


5 


-0.008 


1.00 


1.00 


349 


5 


-0.008 


0.67 


0.45 


355 


5 


-0.011 


0.65 


0.34 


360 


5 


0.002 


0.58 


0.31 


375 


5 


0.015 


0.53 


0.24 


341 


1 


0.007 


1.00 


1.00 


349 


1 


0.004 


0.60 


0.08 


355 


1 


0.0004 


0.60 


0.21 


360 


1 


0.0002 


0.53 


0.15 


375 


1 


0.005 


0.50 


0.16 



the variations in polarized intensity and polarization angle 
must be due to a combination of two processes: Faraday 
rotation and depolarization. We will return to these in 
Sects. H and 
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Fig. 3. Five maps of a part of the Auriga region at frequencies 341, 349, 355, 360 and 375 MHz. The grey scale 
is polarized intensity P, oversampled 5 times, where white is high intensity. Maximum P (in the 375 MHz band) is 
approximately 13 K. Superimposed are polarization pseudo- vectors, where each line denotes an independent synthesized 
beam. 



3.3. Polarization angle 

In Fig. |21 we show, for a 3°x2° subfield slightly northwest 
of the center of the mosaic, maps of polarized intensity 
P (grey scale) and polarization angle <f> (superimposed 
pseudo- vectors) in the 5 frequency bands. In the grey 
scale, which has 5 samples per linear beam width, white 
corresponds to high P. The length of the vectors (one 



per beam) scales with polarized intensity. The amount of 
structure in polarization angle is highly variable. At high 
polarized intensities, the polarization angles mostly vary 
quite smoothly. However, there are also abrupt changes on 
the scale of a beam, the most conspicuous of which give 
rise to the depolarization canals. 
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Fig. 4. Two typical parts of 3x5 beams in the observed field, represented in small graphs of polarization angle <f> 
against A 2 , one graph per independent synthesized beam. The left 3x5 graphs are taken from a part of the field with 
high polarized intensity (<j>15 per frequency band), the right part from a region of lower intensity [a « 5). The 
values within parentheses above each graph are RM and reduced % 2 , respectively. In general, RMs are more regular 
and alike over several beams in high polarized intensity regions, but there is coherence in RM in lower P regions as 
well. 



4. Faraday rotation 

4.1. The derivation of rotation measures 

The rotation measure (RM) of the Faraday-rotating mate- 
rial follows from 4>(\ 2 ) oc RM A 2 . As <j> = 0.5 arctan(/7/<2) 
2 , the values of polarization angle are ambiguous over 
±n 180° (n = 1,2,3,.. .) because a value of angle(Q, U) 
is ambiguous over ±n 360°. For most beams, it was possi- 
ble to obtain a linear c/>(A 2 )-relation by using angles with a 
minimum angle difference between neighboring frequency 
bands. Only in about 1% of the data, addition or sub- 
traction 180° to the polarization angle in at least one 
of the frequency bands was needed. One could increase 
RMs by adding or subtracting multiples of 180°, but that 
would lead to RMs in excess of 100 rad m~ 2 . If such high 
RMs were real, bandwidth depolarization (see Sect. 0) 
would have totally annulled any polarized signal. In addi- 
tion, the RMs we obtain using the criterion of minimum 
angle differences (of order of 10 rad m -2 ) were also ob- 
tained in other observations in this direction (Bingham & 
Shakeshaft 1967, Spoelstra 1984). So we conclude there 
are hardly any 180°-ambiguities in polarization angle and 



2 See footnote in Sect.|2] 



we use the minimum angle difference between frequency 
bands to determine RMs. 

In Fig. 0] we show </>(A 2 )-plots for small arrays of con- 
tiguous, independent beams in typical subfields of the 5'- 
resolution map. In parentheses are the values of the RM 
and the reduced \ 2 °f the linear 0(A 2 )-fit. The leftmost 
3x5 plots are from a region with high P (<r> 15 in all 
frequency bands), while the rightmost plots are from a re- 
gion with lower P (a « 5). In general, the 0(A 2 )-relation 
is closer to linear in high P regions than in lower P re- 
gions, but in the larger part of the Auriga field RMs are 
very well-determined. 

4.2. The distribution of rotation measures 

The RMs in the Auriga- field are plotted as circles in 
Fig. El Filled circles denote positive values, open circles 
negative, and the diameters of the circles are proportional 
to RM. For clarity, we only show the RMs for one in 
four independent beams. In addition, RMs are only shown 
when they were reliably determined, i.e. for the average 
polarized intensity (P) > 5<t, and for reduced \ 2 of the fit 
< 2. 

The left panel shows RM as determined from the ob- 
servations. Most RMs are negative; the few positive values 
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Fig. 6. Histogram of RMs in the Auriga field. The dashed 
line gives the distribution of all "reliably determined" 
RMs, i.e. where P > 5a and reduced x 2 < 2. The solid 
line gives the reliably determined RM distribution after 
subtracting a best-fit linear RM gradient from the data. 
Dotted lines show Gaussian fits to the two histograms. 

occur for 5 ^ 52°, while the largest negative values occur 
for 5>55°. We modeled this systematic variation by a 
large-scale linear gradient, and find a gradient of about 
1 rad m -2 per degree, with the steepest slope along posi- 
tion angle (from north through east) of —20°. 

The average RM, i.e. the RM value in the center of the 
fitted i?M-plane, over the field is RM « -3.4 rad m~ 2 . 
Subtracting the best-fit gradient from the RM distribu- 
tion yields the RMs shown in the right-hand panel in 
Fig. with only small-sale structure in RM. The distri- 
bution of small-scale RMs (i.e. RM where the large-scale 
gradient is subtracted) is significantly narrower than the 
total RM distribution (<7 RM = 1.4 K instead of 2.3 K), and 
is more symmetrical, as is shown in Fig. |SJ where the two 
histograms are compared. 

The decomposition of the i?Af-distribution into a con- 
stant component, a gradient, and a small-scale component 
is not physical, as there is probably structure on all scales. 
However, for our purpose it is a good approximation, and 
it allows us to estimate the large-scale and random com- 
ponents of the Galactic magnetic field, see Sect. [7] 

4.3. Structure in RM on arcminute scales 

Although the RM -distribution shows structure on many 
scales, the most intriguing changes occur on the scale of 
the beam, as illustrated in Fig.0 In the figure, we display 
a 9x7 array of </>(A 2 )-plots, overlaid on a grey scale repre- 
sentation of P at 349 MHz. Although RM varies quite 
smoothly over most of the area, there are also abrupt 
changes, which frequently involve a change of the sign of 
RM (e.g. at (a, 8) = (94.60, 53.00)° or (94.68, 53.20)°). 
Abrupt RM changes over one beam with the right mag- 
nitude to cause an angle change of ±(n+ 1/2) 180° cause 
depolarization canals. As RM is an integral over the en- 
tire line of sight, it is difficult to understand such changes 
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Fig. 7. White symbols and lines are plots of 4> against 
A 2 and their linear fits, one for each independent beam, so 
that the slope is RM. Sudden RM changes occur over one 
beam width. The grey scale is polarized intensity at 349 
MHz oversampled by a factor of 5. Maximum P (lightest 
grey) is ~ 40 mJy/beam, while in the bottom left pixel P 
is set to zero for comparison. 

in RM as resulting from structure in the magnetic field 
and/or distribution of electrons. As a matter of fact, nu- 
merical simulations of the warm ISM show that beam 
depolarization can enhance RM differences and greatly 
steepen the gradient (Haverkorn 2002). In addition, mod- 
els of a medium that produces Faraday rotation and that 
emits synchrotron radiation show that the RM can change 
sign without a corresponding change in the direction of the 
magnetic field (Sokoloff et al. 1998, Chi et al. 1997). 

5. The structure in P, and the implied properties 
of the ISM 

As / and P are uncorrelated, the structure in P cannot 
be caused by structure in emission, but instead is created 
by several instrumental and/or depolarization effects. We 
discuss these processes in detail in two forthcoming papers 
for several fields of observation, among which the Auriga 
field (Haverkorn et al. 2003b, c). Here we summarize them, 
and discuss the results for the Auriga region. 

Structure in P could be induced by the insensitivity 
of the WSRT to large-scale structure, as a result of miss- 
ing short-spacing information. However, it is very unlikely 
that this effect is important in the Auriga observations, 
because the distribution of RMs per pointing center is 
too wide to allow offsets. If the distribution of RMs in a 
Faraday screen is too wide, the distributions of Q and U 
of a polarized background that is Faraday-rotated by the 
screen have averages very close to zero (even if the back- 
ground was uniformly polarized). As a result, the missing 
short-spacing information does not produce offsets. For 
our frequencies the latter is true if crm>1-8 rad m~ 2 . 
In the Auriga region, this condition is met in essentially 
all of the subfields of the mosaic. Furthermore, a test in 
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Fig. 5. Rotation measures given as circles with a diameter proportional to the magnitude of the RM , superimposed 
over P at 349 MHz in grey scale. The scaling is in rad m~ 2 . RMs are only shown if P > 5<r and reduced \ 2 < 2, 
and in both directions, only every second independent beam is plotted for clarity. The right map shows the same data 
where the best-fitting linear gradient in RM is subtracted. 



which offsets are added (described in the papers referred 
to above) to produce better linearity of the </>(A 2 )-relation, 
did not result in reliable non-zero offsets. 

Structure in P could also be generated by non- 
constant, but high, RM values as a result of bandwidth 
depolarization. Summation of polarization vectors with 
different angles in a frequency band will reduce the polar- 
ized intensity. The amount of depolarization depends on 
the spread in polarization angle within the band which, 
in turn, depends on the RM and the width of the fre- 
quency band. For the RMs of order 10 rad m -2 and our 
bandwidth of 5 MHz the effect is negligible. 

A special kind of structure in P is due to beam depo- 
larization, viz. the cancellation of polarization by vector 
summation across the beam defined by the instrument. 
This is a 2-dimcnsional summation of polarization vectors 
across the plane of the sky, which is probably responsible 
for the depolarization canals. It certainly cannot produce 
all of the structure in P because it occurs on a single, 
small scale, whereas P has structure on a large range of 
scales. 

Finally, depth depolarization occurs along the line of 
sight in a medium that both emits synchrotron radiation 
and produces Faraday rotation. It must be responsible 
for a large part of the structure in P. The ISM in the 
Galactic disk most likely is such a medium. The Galactic 
synchrotron emission has been modeled by Bcucrmann et 
al. (1985) as the sum of two contributions: the thin disk 
with a half equivalent width of ~ 180 pc, which coincides 



with the stellar disk, and a thick disk of half equivalent 
width ~ 1800 pc (where the galactocentric radius of the 
sun is assumed Rq — 10 kpc). Thermal electrons exist 
also out to a scale height of approximately a kpc, con- 
tained in the Reynolds layer (Reynolds 1989, 1991). So up 
to a height of about a kpc, the medium probably consists 
of both emitting and Faraday-rotating material. In such a 
medium, the polarization plane of the radiation emitted at 
different depths is Faraday-rotated by different amounts 
along the line of sight. Therefore the polarization vectors 
cancel partly, causing depolarization along the line of sight 
even if the magnetic field and thermal and relativistic elec- 
tron densities are constant (differential Faraday rotation, 
see Gardner & Whiteoak 1966, Burn 1966, Sokoloff et al. 
1998). If the medium contains a random magnetic field, 
the polarization plane of the synchrotron emission will 
vary along the line of sight causing extra depolarization in 
addition to the differential Faraday rotation (called inter- 
nal Faraday dispersion). We will refer to the combination 
of the two mechanisms along one line of sight as depth de- 
polarization, thus adopting the terminology proposed by 
Landecker et al. (2001). 

Depth depolarization is the dominant depolarization 
process in the Auriga field. We constructed a numerical 
model in Haverkorn et al. (2003c) to evaluate if depth de- 
polarization can indeed be held responsible for most of 
the structure in P, in the Auriga region as well as in the 
other region that we studied. At the same time, the model 
makes it possible to estimate the conditions in the warm 
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component of the ISM. The model consists of a layer cor- 
responding to the thin Galactic disk, where random mag- 
netic fields on small scales are present, and a constant 
background polarized radiation. 

In the model, the thin disk (with height h) is divided 
in cells of a single size d, which contain a random magnetic 
field component B ran and a regular magnetic field B reg . 
The amplitude of B ran is constant, but in each cell its 
direction is randomly chosen. Each cell emits synchrotron 
radiation, proportional to B\ = (B ran ,±_) 2 + {B re g,±) 2 , 
where we add energy densities because the two compo- 
nents are physically separated different fields. Only a frac- 
tion / of the cells contains thermal electrons, and there- 
fore Faraday-rotates all incoming radiation (to represent 
the filling factor / of the warm ISM). 

We distinguish several kinds of model parameters. 
First, we have input parameters with values taken from 
the literature; these are thermal electron density n e = 
0.08 cm~ 3 , filling factor / = 0.2 (both from Reynolds 
1991), and height of the Galactic thin disk h — 180 pc 
(Beuermann et al. 1985). Secondly, there are constraints 
from the observations of the Auriga field, viz. the width, 
shape and mean of the Q, U, I and RM distributions. In 
addition, we define parameters without observational con- 
straints, such as cell size d and emissivity per cell, which 
we vary in a attempt to constrain them from the model 
given the observational constraints. Finally, there are pa- 
rameters that can be adjusted and optimized for each al- 
lowed cell size and emissivity per cell, like the random 
and regular magnetic field components and background 
intensity. 

Using the observations of the Auriga region to con- 
strain the parameters in the depth depolarization model, 
we obtain the following estimates. The cell size d is prob- 
ably in the range from 5-20 pc, the random component 
of the magnetic field B ran is about 1 /xG for a cell size of 
5 pc or larger, and up to 4 /iG for smaller cell sizes. With 
the estimate of the large-scale, regular field (<~ 3 fiG) this 
leads to a value of about 0.3 to 1.3 for B ran / B reg . This 
is lower than most other estimates of B ran / B reg in the 
literature, which seem to indicate a value larger than 1. 
One possible explanation for this discrepancy is the fact 
that the Auriga region is not very large. Therefore ran- 
dom magnetic field components that project to angular 
scales larger than a few degrees will have been interpreted 
as components of the regular field. If such large-scale ran- 
dom components exist, this effect will artificially decrease 
Bran/ B reg . Secondly, the Auriga field was selected for its 
conspicuous polarization structure, and high polarized in- 
tensity in general indicates a more regular magnetic field 
structure. Finally, the Auriga region is observed in an 
inter-arm region. The next spiral arm, the Perseus arm, 
is more than 2 kpc away, and we expect that at these fre- 
quencies, all polarized emission from the Perseus arm is 
depolarized. There are indications that in the inter-arm re- 
gion, the regular magnetic field component is higher than 
the random magnetic field component (Han & Qiao 1994, 
Indrani and Deshpande 1998, Beck 2001). 
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Fig. 8. The amount of polarization contributed by a cer- 
tain cell, against the distance to that cell, in a single-cell- 
size model of polarized radiation propagating through a 
magneto-ionized medium. 

A cell of size 10 pc located at the far edge of the thin 
disk still subtends a little less than a degree in our maps. 
We see structure on smaller scales, so that cells of approx- 
imately a parsec have to be present as well. Indeed, there 
are many indications that structure exists in the ISM on 
many scales (Armstrong et al. 1995, Minter & Spangler 
1996). 

The distance out to which polarized radiation is not 
depolarized by depth depolarization, as derived from the 
model, is shown in Fig. |S] The total path length through 
the medium containing small-scale structure was about 
600 pc, which corresponds to a height of the thin disk of 
about 150 pc for the Auriga latitude of ~ 16°. Radiation 
emitted at large distances is depolarized more than emis- 
sion from nearby, but there is no definite polarization hori- 
zon due to depth depolarization. Beam depolarization can 
cause a polarization horizon, as structure further away will 
be on smaller angular scales, which causes more beam de- 
polarization (Landecker et al. 2001). 

6. Polarized extragalactic point sources 

Thirteen polarized extragalactic point sources in the 
Auriga-field were selected based on three selection crite- 
ria: (1) signal-to-noise a > 3 in P at all frequencies, (2) 
the degree of polarization p is higher than 1% at all fre- 
quencies (because a lower polarization degree can also be 
caused by instrumental polarization, see Sect.|2J), and (3) 
the source does not lie too far at the edge of the field away 
from any pointing center, to limit the contribution from 
instrumental polarization. The 13 sources are all detected 
in the 1' resolution data, and most of them are unresolved 
at that resolution. Their RMs vary from approximately 
-13 to 13 rad m~ 2 . 

The relevant properties of the sources are given in 
Table |3 where the second column gives the right ascen- 
sion and declination, column 3 the RM, and column 4 
the reduced \ 2 °f the </)(A 2 )-fits. Columns 5, 6, and 7 give 
the values of P, I and percentage of polarization p respec- 
tively, averaged over all frequency bands. Fig. shows 
the polarization angle <p plotted against A 2 for each point 
source. In Fig. we show the 13 sources overlaid on a 
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Table 3. Extragalactic sources in the Auriga region. 
The second column gives position, and the third RMs. 
Reduced x 2 of the </>(A 2 ) relation is given in column 4. 
Columns 5, 6 and 7 give resp. P, / (both in mJy/beam) 
and percentage of polarization p averaged over the five 
frequency bands. 



No. 


°) 


RM (rad m" 2 ) 


x 2 


(P) 


(I) 


(P) 


1 


[ 97.9, 55.3 ] 


7.4 ±0.2 


12.9 


45.6 


2021 


2.3 


2 


[ 97.7, 53.6 ] 


12.9 ± 0.3 


3.1 


14.9 


216 


6.9 


3 


[ 96.6, 53.3 ] 


6.6 ±0.4 


1.7 


15.0 


224 


6.7 


4 


[ 96.5, 54.1 ] 


5.7 ±0.4 


3.7 


11.9 


380 


3.1 


5 


[ 95.3, 50.8 ] 


0.6 ±0.4 


1.6 


14.7 


165 


8.8 


6 


[ 95.2, 50.8 ] 


3.1 ±0.4 


0.2 


14.7 


197 


7.5 


7 


[ 94.7, 54.8 ] 


7.9 ±0.3 


3.0 


18.8 


282 


6.7 


8 


[ 93.4, 55.6 ] 


6.3 ±0.2 


1.1 


21.3 


180 


11.8 


9 


[ 92.3, 52.9 ] 


2.7 ±0.1 


6.1 


35.8 


1222 


2.9 


10 


[ 90.3, 50.4 ] 


-10.2 ± 1.4 


4.4 


7.4 


168 


4.4 


11 


[ 89.3, 55.0 ] 


3.0 ±0.4 


1.4 


11.6 


340 


3.4 


12 


[ 89.2, 52.5 ] 


-8.1 ±0.1 


12.7 


22.1 


477 


4.6 


13 


[ 89.0, 52.3 ] 


-13.6 ±0.4 


5.8 


14 


237 


5.9 



grey scale plot of polarized intensity at 349 MHz, with 
their RMs indicated. The diameter of the symbol and its 
shading indicate the value of RM. 

Sources 1 and 2 could have a significant contribution 
of instrumental polarization, because they are observed in 
only one pointing center and are located about 0.8° to 1° 
away from the pointing center. Nevertheless, they show 
linear </>(A 2 )-relations, which indicates that instrumental 
polarization is not important. These sources are therefore 
included in the analysis. 

The strong correlation of RM across the sky indicates 
a Galactic component to the RMs of the extragalactic 
point sources. The best-fitting linear gradient to the RMs 
of the sources has a steepest slope of 3.62 rad m -2 per 
degree in position angle 72°, i.e. roughly in the direc- 
tion of increasing Galactic latitude (position angle 66°). 
The standard deviation of RMs around this gradient is 
ctrm ~ 3.6 rad m -2 . This is consistent with the result of 
Leahy (1987), who finds a typical 'internal' RM contribu- 
tion from the extragalactic source or a halo around it of 
~ 5 rad m~ 2 . 

The change of sign of RM of the extragalactic sources 
within the observed region indicates a reversal in the 
Galactic magnetic field. (Unlike the diffuse emission, sign 
changes in RMs of extragalactic sources cannot be due 
to depolarization effects.) The reversal exists on scales 
larger than our field of view, as is shown in Fig. 1111 where 
we combine the RMs of our sources with those from the 
literature. The circles are our sources, the squares indi- 
cate sources that were detected by Simard-Normandin et 
al. (1981) and/or by Tabara & Inoue (1980), and the 
one triangle indicates the only pulsar nearby (Hamilton 
& Lyne 1987). The numbers next to the squares and 
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Fig. 9. Plots of polarization angle 4> against A 2 for the 13 
polarized extragalactic sources in Table El Note that the 
scaling of the (/>-axis differs for each source. 
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Fig. 10. RMs of the 13 polarized extragalactic point 
sources in the Auriga field in circles, overlaid on P at 
349 MHz. Filled circles are positive RMs, open circles 
negative. The crosses denote the pointing centers, and the 
scaling is in rad m . 
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Fig. 11. RMs of extragalactic point sources and of one 
pulsar. The circles are the point sources detected in the 
Auriga field, where the radius of the circle is propor- 
tional to the magnitude of RM. Squares are extragalactic 
point sources detected by Simard-Normandin et al. (1981) 
and/or by Tabara & Inoue (1980), and the one triangle 
in the field denotes the only pulsar near the Auriga field 
(Hamilton & Lyne 1987). Values of the RM are written 
next to the sources, and for RM < 15 rad m -2 the radius 
of the symbol is proportional to RM. 



triangle give the magnitude of RM of the source (the 
diameter of the symbol is only proportional to RM 
up to RM = 15 rad vcT 2 ). One source at (a, 6) = 
(91.25°, 48°) was omitted, because Simard-Normandin et 
al. give RM = 34 rad m -2 and Tabara & Inoue give 
RM = —64.9 rad m -2 . The source distribution shows 
a clear magnetic field reversal, which is not on Galactic 
scale but on smaller scales, as is shown in Figs. 1 and 2 of 
Simard-Normandin & Kronberg (1980), which show RMs 
of extragalactic sources over the whole sky. 

The gradient in RM measured in the extragalactic 
point sources is completely unrelated and almost perpen- 
dicular to the structure in RM from the diffuse emis- 
sion. This can be explained from the widely different path 
lengths that extragalactic sources and diffuse emission 
probe. Diffuse emission can only be observed out to a dis- 
tance of a few hundred pc to a kpc, as radiation from fur- 
ther away will be mostly depolarized. On the contrary, ex- 
tragalactic sources are Faraday-rotated over the complete 
path length through the Milky Way of many kiloparsecs 
long. Therefore, the RM structure from diffuse radiation 
and from extragalactic sources gives information about 
different regions in the Galaxy. 



7. Discussion 

7.1. High P -structures and alignment with the 
Galactic plane 

The depth depolarization model described in Sect. only 
describes the properties of the ISM in a statistical man- 
ner, by giving a representation of an average line of sight. 
Therefore, it does not pretend to give a detailed descrip- 
tion of the distribution of P, as it is actually observed. 
Even the global properties of the P distribution, like the 
general alignment with the Galactic plane, is not part 
of the model. This alignment is not only visible in the 
Auriga field, but also in other regions observed with the 
WSRT (Haverkorn et al. 2003a, Schnitzeler et al. in prep) . 
The large-scale structures in P cannot be caused by sub- 
stantial extra emission, because that should be accompa- 
nied by corresponding structure in I, which is not ob- 
served. Instead, where the polarized intensity is highest, 
the depolarization is probably the least, so that there the 
ISM is relatively transparent to polarized emission. In 
general, high depolarization is caused by a large amount 
of structure in RM, along the line of sight and/or over 
the sky (the latter on the scale of the beam or smaller). 
This explanation of the structure of high P is borne out 
by the relation between polarized intensity P and er RM , 
shown in Fig. [El This fi gure shows the width (T RM of a 
Gaussian fit to the RM distribution in several intervals 
of polarized intensity with a width AP — 0.01 Jy/beam. 
Rotation measure clearly varies more at lower polarized 
intensity, where only well-determined RMs are used so 
that the effect cannot be due to noise. For the two in- 
tervals of highest P (0.08 - 0.09 Jy/beam and 0.09 - 
0.10 Jy/beam), not shown in Fig.^1 no Gaussian fit could 
be made, as the distributions were bimodal. The peaks in 
these bimodal distributions correspond to two spatially 
coherent structures, which are shown in the two panels 
of Fig. CHI In the region around (a, 6) = (94.6°, 51.8°), 
— 1.5 Ss RM ^ rad m -2 , while in the region around 
(a, S) = (92.2°, 53.8°), <, RM <, 1.3 rad m~ 2 . Clearly, 
RM is indeed very uniform in regions of high polarization. 
The low variation in RM is also reflected in the uniformity 
of the polarization angles. 

A constant RM over a certain area sets an upper limit 
on structure in magnetic field and thermal electron den- 
sity. The fact that some of the filamentary structure in 
P is aligned with Galactic latitude suggests a magnetic 
origin, although stratification of n e over Galactic latitude 
can contribute as well. For variation of polarization angle 
over the filament of less than about 20° (from Fig. I13|) , 
ARM ^ 0.5 rad m -2 is needed. If magnetic field struc- 
ture would be the same over the entire path length towards 
the filament, and assuming n e = 0.08 cm -3 , this requires 
a magnetic field change ABu <> 0.01 /xG with a polariza- 
tion depth of 600 pc. Even if the filamentary structures 
are sheet-like with a constant magnetic field over a large 
part of the path length, or if the polarization horizon is 
locally much shorter than the highly uncertain value of 
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Fig. 12. Width of a Gaussian fit to the RM distribu- 
tion (T RM for a range of P, where P is divided in bins 
of 0.01 Jy/beam. 
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Fig. 13. Two areas in the P map at 349 MHz with high- 
est P (> 80 m Jy/beam). Overlaid pseudo- vectors are po- 
larization angles, showing a minimal change in angle for 
high P. 

600 pc, the upper limit to AB is very low. Therefore, 
this structure of high P aligned with Galactic latitude 
indicates the existence of long filaments or sheets paral- 
lel to the Galactic plane of highly uniform magnetic field 
and thermal electron density. As back- and/or foreground 
structure can also imprint structure in RM , these regions 
may well be characterized by low thermal electron density. 

7.2. The Galactic magnetic field 

The Auriga field of 7°x9° represents only 0.1% of the 
sky, and is therefore by itself not very well suited for an 
analysis of the global properties of the Galactic magnetic 
field. However, we can estimate strength and direction of 
the regular component of the Galactic magnetic field from 
the distribution of RM s of the diffuse radiation and of the 
extragalactic sources over the field. 

First, we need to determine how much of the RM gra- 
dient seen in the diffuse emission is caused by n e , and how 
much by B». Independent measurements of n e can be ob- 
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Fig. 14. Superimposed circles are emission measures from 
the WHAM Ha survey, to a maximum intensity of 
~ 5 Rayleigh, at a resolution of one degree. The grey scale 
denotes P at 349 MHz. 



tained from the emission measure EM = J n 2 e dl, as mea- 
sured with the Wisconsin Ha Mapper survey (WHAM, 
Haffner et al., in prep, Reynolds et al. 1998), shown in 
Fig. O As the resolution is about a degree, the WHAM 
survey cannot be used to determine small-scale structure 
in n e , but gives an estimate of the global gradient in n e 
over the field. From north to south, the Ha observations 
show an increase from about 1 to 5 Rayleighs (1 Rayleigh 
is equal to a brightness of 10 6 /47r photons cm~ 2 s _1 
ster -1 , and corresponds to an EM of about 2 cm -6 pc for 
gas with a temperature T = 10000 K). The small enhance- 
ment in EM at (a,S) = (94.8°, 55.5°) is probably related 
to the ancient planetary nebula PuWel (PW1, Tweedy & 
Kwitter 1996). Note that the Ha emission probes a longer 
line of sight than the diffuse polarization, and is velocity 
integrated over —80 to 80 km s _1 LSR, so complete cor- 
respondence is not expected. However, a gradient in RM 
that was solely due to the n e -gradient would have a sign 
opposite to that of the observed i?M-gradient. Therefore, 
the observed gradient in RM must be due to a gradi- 
ent in the parallel component of the magnetic field. If the 
gradient in Ha originates in the same medium as the po- 
larized radiation, the gradient in magnetic field must even 
be large enough to compensate the gradient in electron 
density in the opposite direction. 

So the RM gradient of the diffuse polarized emission 
of about 1 rad m~ 2 per degree in the direction of posi- 
tive Galactic longitude must be caused by a corresponding 
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increase of the parallel component of the Galactic mag- 
netic field. The direction of the gradient in longitude and 
the sign of the average RM are consistent with a regu- 
lar Galactic magnetic field directed away from us in the 
second quadrant, as found from other observations (see 
e.g. review by Han, 2001). Attributing the observed RM- 
gradient completely to the change in the parallel com- 
ponent of the Galactic magnetic field, we can derive the 
strength of its regular component. The average RM$ and 
the value of the RM gradient give two independent esti- 
mates of the regular magnetic field, if a certain pitch angle 
is assumed. Assuming a path length of 600 pc and electron 
density n e — 0.08 cm~ 3 , the two magnetic field determi- 
nations agree on a regular magnetic field of B reg « 1 fiG 
for a pitch angle p = —14°, in agreement with earlier es- 
timates (Vallee 1995). 

The gradient in the RMs of the extragalactic sources 
does not show a longitudinal component at all, contrary 
to that in the RMs of the diffuse polarization observa- 
tions. This indicates that other large-scale magnetic fields 
become important, fields that vary on scales larger than 
our field. Observations of halos of external galaxies (e.g., 
Dumkc ct al. 1995) have shown that cell sizes in the 
Galactic halo are much larger than in the disk, from 100 
- 1000 pc. Even at a distance of 3000 pc, the size of the 
total Auriga field is ~ 400 pc, so big cells in the halo could 
be responsible for a considerable part of what we call the 
regular component of the magnetic field, as measured in 
the RMs of the extragalactic sources. 

We subtracted the gradient of the diffuse emission from 
the gradient in the RMs of the extragalactic sources, af- 
ter which a gradient in RM of ARM « 3.6 rad m~ 2 in 
position angle 54° remains (only deviating by about 10° 
from the direction perpendicular to the Galactic plane). 
This indicates a regular magnetic field of about 1 /iG in 
this direction, assuming that n e and B remain constant 
over a total line of sight of 3000 pc. The magnetic field 
strengths actually present in the medium will be higher if 
the field varies along the line of sight. 

8. Conclusions 

Multi-frequency polarization observations of the diffuse 
Galactic background yield information on structure in the 
Galactic warm gas and magnetic field. 

The multi-frequency WSRT observations of a region 
of size ~ 7°x9° at I = 161°, b = 16° in the constellation 
Auriga show a smooth total intensity /, but abundant 
structure in Stokes parameters Q and U on several scales, 
with a maximum Tb_ po i ~ 13 K (about 30% polarization). 
Filamentary structure up to many degrees long is present 
in P, sometimes aligned with Galactic latitude. In ad- 
dition, narrow, one-beam wide depolarization canals are 
most likely created by beam depolarization and may well 
indicate abrupt RM changes. As the polarization struc- 
ture is uncorrelated to /, the structure in P cannot be cre- 
ated by variations in synchrotron emission, but has to be 
due to Faraday rotation and depolarization mechanisms. 



Rotation measure maps show abundant structure on 
many scales, including a linear gradient of <~ 1 rad m~ 2 
per degree in the direction of positive Galactic longitude 
and an average RM ~ —3.4 rad mT 2 . The gradient is 
consistent with the regular Galactic magnetic field if its 
strength is B reg « 1 [iG and the field is azimuthally ori- 
ented with a pitch angle p « —14°. Ubiquitous structure is 
present in the RM map on beam size scales (5'), indicat- 
ing significant changes in magnetic field and/or thermal 
electron density over very small spatial scales. 

There are two dominant depolarization mechanisms 
which create structure in P in the Auriga field. First, 
beam depolarization (i.e. depolarization due to averaging 
out polarization angle structure within one beam width) 
most likely creates the canals, and is important in re- 
gions of low P. Additional depolarization is needed to ex- 
plain the observed P distribution, which is only possible if 
the medium both Faraday-rotates and emits synchrotron 
emission. Then, emission originating in the medium can be 
depolarized along its path towards the observer, so-called 
depth depolarization. The polarization angle along the 
path can vary due to Faraday rotation, or due to varying 
intrinsic polarization angle of the emission. This indicates 
a varying magnetic field, which however is constrained by 
the upper limit on structure in /. A depth depolariza- 
tion model was constructed of a layer of cells with varying 
magnetic fields and a constant background polarization. 
Constraints from the Auriga observations give estimates 
for several parameters in the ISM. The cell size of struc- 
ture in the ISM is constrained to ~ 15 pc, and the ratio 
of random to regular field is 0.7 ± 0.5. 

Thirteen extragalactic sources in the Auriga field also 
show a gradient in RM, but roughly in the direction of 
Galactic latitude, which is perpendicular to the gradient 
in the diffuse emission. The RM distributions from dif- 
fuse radiation and from extragalactic point sources are so 
different because the diffuse radiation mostly probes the 
first few hundred parsecs, whereas RMs from the point 
sources are built up along the entire line of sight through 
the Milky Way. RMs of extragalactic sources change sign 
over the field, which indicates a local reversal of the mag- 
netic field. 
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